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abstract 


BLAST  LOADING  OF  MULTI-STORIED  STRUCTURES 
OF  FIRMO- VISCOUS  TYPE 

by 

W.H.  Hoppmann  II 
N.J.  Huffington,  Jr. 


This  study  is  the  second  in  a general  analytical  investigation  of 
the  blast  loading  of  structures.  Previously,  the  various  elements 
of  the  problem  were  enumerated  and  discussed.  Here,  analyses  are 
presented  of  the  deformation  of  two  types  of  structures  assumed  to 
have  internal  resistance  of  the  firtno-viscous  type  and  to  be  sub- 
jected to  side  blast  loading.  The  first  is  a multi-storied  building 
consisting  of  a series  of  massive  floors  each  of  which  is  supported 
atop  shear  frames  connecting  the  successive  floors.  The  second 
is  a pile  of  rigid  blocks  cemented  togethor,  intended  to  simulate 
massive  blocks  of  masonry  attached  to  each  other  by  a viscous 
material.  Numerical  examples  illustrating  the  theory  are  treated. 
This  study  indicated  that: 

The  type  of  analysis  presented  may  be  of  value  in 
• predicting  damage  caused  by  blast  against  a multi- 
storied building  or  similar  structure. 

^ There  is  a need  for  more  reliable  information  on 
structural  damping. 
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SUMMARY 


PROBLEM 

To  analyze  dynamically  the  deformation  of  two  types  of  structures  which  are 
assumed  to  have  internal  resistance  of  the  firmo-viscous  type  and  to  l>e  subjected  to  side 
blast  loading. 

FACTS 

This  study  is  the  second  in  a general  analytical  investigation  of  the  blast  loading 
of  structures.  Previously,  the  various  elements  of  the  problem  were  enumerated  and 
discussed.  Among  these  were  stress-strain  laws,  including  damping.  An  illustrative  prob- 
lem of  internal  forces  was  analyzed  The  dynamical  analysis,  however  defective,  is  supe- 
rior to  surmises  of  structural  damage  based  on  statical  analysis  only  and  is  amenable  to 
simple  numerical  calculation 

DISCUSSION 

In  various  tall  buildings  the  resistance  to  applied  loads  may  be  essentially  of  a 
shear  type  Also,  there  may  be  a viscous  element  in  the  total  internal  resistance.  A method 
for  analyzing  these  cases  is  presented  A structure  is  supposed  subjected  to  side  blast  of 
the  type  described  in  “The  Effects  of  Atomic  Weapons”  (3).  This  structure  is  representable 
as  a series  of  rigid  masses  supported  by  columns  which  are  not  cross-braced.  Another 
structural  model  of  considerable  interest  is  that  of  a vertical  system  of  rigid  masses  acted 
upon  by  external  impulsive  forces  and  coupling  forces  which  are  internal  to  the  structure. 
This  model  is  especially  applicable  to  structures  composed  of  massive  blocks  of  masonry 
cemented  together;  it  also  has  application  to  buildings  for  which  the  mass  may  reasonably 
be  assumed  to  be  concentrated  at  floor  levels.  The  latter  application  was  suggested  in 
“The  Effects- of  Atomic  Weapons”  (3)  and  was  used  earlier  for  a building  considered  as 
a two-degrce-of-freedom  system  subjected  to  earthquake  motions. 


CONCLUSIONS 

I 1 he  type  of  analysis  presented  is  considered  to  be  of  value  in  predicting  damage 
caused  by  blast  against  a multi-storied  building  or  similar  structure.  This  point  of  view 
is  advanced  despite  obvious  defects  in  idealizations  of  the  type  involved. 

2.  The  method  of  treating  shear-frame  building  may  be  utilized  in  the  analysis  of 
similar  structural  units  occurring  in  equipment  and  machines  subjected  to  dynamic  loads 
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However,  it  will  still  lx'  necessary  to  determine  suitable  stress-strain  laws  in  each  particular 
case  of  interest. 

3.  Investigations  involving  the  firmo-viscous  hypothesis  emphasize  the  need  for 
more  reliable  information  on  structural  damping.  Considerable  effort  to  study  the  use  of 
damping  laws  in  structural  analysis  seems  to  be  warranted  at  the  present  time. 
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INTRODUCTION 


In  a previous  paper  (2)  the  various  elements  of  the  problem  of  blast  loading  of 
structures  were  enumerated  and  discussed.  Among  these  were  stress-strain  laws,  including 
damping.  Also,  in  connection  with  an  illustrative  problem  which  was  solved  it  was  noted 
that  the  formal  mathematical  description  and  solution  might  represent  a case  in  which 
resisting  forces  were  of  the  “perfectly  plastic”  type  or  equally  well  a case  in  which  the 
resisting  forces  were  of  the  Coulomb  frictional  type.  While  the  mathematical  equations 
describing  these  two  cases  are  identical  the  physical  reality  considered  is  different.  In  the 
one  the  resisting  force  is  produced  by  the  plastic  deformation  of  a coupling  material  be- 
tween the  successive  blocks  composing  the  idealized  structure.  In  the  other  case  the  re- 
sisting force  is  thought  of  as  produced  by  the  rubbing  together  of  two  adjacent  surfaces. 
The  first  case  is  a modification  of  that  treated  in  “The  Kffects  of  Atomic  Weapons”  (3). 

Many  other  kinds  of  internal  resistance  are  treated  in  rheology  (1)  which  is  the 
science  of  the  deformation  and  flow  of  matter.  Some  of  the  laws  of  resistance  are  shown 
in  Fig.  1.  They  are  represented  on  three  axes  representing  the  shearing  stress  r,  the  shear- 
ing strain  7,  and  the  time  rate  of  change  of  shearing  strain  7,  respectively.  The  equations 
corresponding  to  the  diagrams  are  also  given.  The  resistance  law  shown  in  Fig.  1(c)  is 
the  one  utilized  in  the  previously  mentioned  paper  (2).  The  one  corresponding  to  the 
Voigt  (5)  or  Kelvin  (4)  solid  is  used  in  the  present  paper.  This  model  has  been  used  con- 
siderably in  structural  analysis  in  the  past  ((>).  It  may  be  considered  to  describe  a damping 
force  in  the  case  of  a vibrating  body  or  simply  as  a distortional  rate  resistance  if  there  is  no 
noticeable  restoring  force  While  the  law  indicates  too  severe  damping  in  the  higher  modes 
of  vibration  of  a structure  (ti)  it  still  may  give  useful  results  in  the  analysis  of  some  vi- 
brating concrete  structures.  Also,  it  may  be  useful  in  predicting  damage  in  structures 
which  are  subjected  to  forces  so  large  that  the  elastic  restoring  forces  may  be  neglected. 
The  general  dynamical  analysis  in  this  latter  case  is  amenable  to  reasonably  simple  nu- 
merical calculation.  This  dynamical  analysis,  however  defective,  is  surely  superior  to  sur- 
mises of  structural  damage  based  on  statical  analysis  only. 

In  various  tall  buildings  the  resistance  to  applied  loads  may  be  essentially  of  a shear 
type.  Also,  there  may  l>e  a viscous  element  in  the  total  internal  resistance.  A method  for 
analyzing  these  cases  is  presented. 
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SIDE  BLAST  LOADING  OF  A MULTI-STORIED  STRUCTURE 

A structure  is  supposed  subjected  to  side  blast  of  the  type  described  in  “The  Effects 
of  Atomic  Weapons”  (3).  This  structure  is  representable  as  a series  of  rigid  masses  sup- 
ported by  columns  which  are  not  cross-braced  (See  Fig.  2.) 


Deformations  are  of  the  shear  rather  than  of  the  bending  type.  This  structure  is 
thought  of  as  a uniform  solid  with  distributed  mass  and  shear  resistance  as  suggested  by 
M.  Biot  (1).  The  mass  density,  effective  section  area,  and  shear  modulus  are  averaged 
quantities.  An  element  of  such  an  idealized  solid  is  shown  in  deformed  position  in  Fig.  3. 
The  shear  frames  deform  so  that 


7 = 


dw 

d~x' 


where  w = displacement  in  transverse  direction. 

Biot  assumed  the  stress-strain  law  to  be  that  of  the  Hooke  solid: 


r = (1 7, 

where  the  shear  modulus  G is  assumed  to  be  an  effective  constant  throughout  the  length. 
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On  this  basis  the  equation  of  motion  can  be  written  in  the  usual  fashion  for  a one-dimensional 
shear  wave : 

G d-w  d-w  . v 

pa7-"¥' 

where  p = average  mass  per  unit  volume. 

The  Biot  derivation  was  intended  for  cases  in  which  the  excitation  of  the  structure 
was  caused  by  earthquake  movement  of  the  ground.  Now  if  a blast  load  P(x,l),  as  shown 
in  Figs.  2-3,  acts  on  the  structure  and  the  stress-strain  law  is  that  of  the  firmo-viscous 
solid  (7),  Fig.  1(d),  expressed  analytically  as 

t = Gy  -I-  P7, 

then  the  equation  of  motion  readily  becomes 

G d-w  p dlw  P(x,t)  _ d*w 

P dx-  P dx 1 dt  pA  dO 

where  p = coefficient  of  viscous  resistance  (assumed  constant) 

A = effective  section  area  in  shear  (assumed  constant). 


Fig.  3 

The  desired  solution  of  the  blast  loading  problem  will  be  a displacement  function 
w(x,t)  which  satisfies  Eq.  1,  the  boundary  conditions 


w = 0 at  x = 0 

AG 

dw  d'w  . 

a,  + ^ a, a, - 0 at 

(2) 

(wjiere  l is  height  of  building), 

and  the  initial  conditions 

w(x,0 ) = 0,  w(x,0)  - 0 

(3) 

for  a structure  which  is  at  rest  before  the  biast  ioading. 
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By  separation  of  variables  and  application  of  the  boundary  conditions  of  Kq.  2 
the  transcendental  frequency  equation  is  found  to  Ik* 


pi  n 
cos  = (), 
a 


where  a - • 


This  equation  is  satisfied  bv  the  roots 

P«  = 


nir  a 
~2l  ’ 


where  n = 1,  3,  5,  . . .. 

Then  the  solution  of  Eq.  1 will  have  the  form 

■£>  nrx 

H = Z,  7-(0  • sln  ~2l  ’ 

n m \ ,3.5  *" 


(4) 


in  which  the  functions  qH(l)  may  be  found  by  use  of  the  Lagrangian  method  (8).  This 
involves  computing  the  kinetic  energy  7’.  the  potential  energy  V,  the  Rayleigh  dissipation 
function  F (9),  and  the  generalized  force  Qn,  which  have  the  values: 

7 = — I u-  <h  = — — 2L 
- •'o  4 B~35 

..  AG  r/dwY-  t t-AG  <-  , 2 

1 - 2 J.(tJ  ‘"'-nr 

Q-  = f^f0  p(j>  sin 

assuming  F(r,t)  = p(r)  ■ f(t). 


Substituting  the  values  of  T,  V,  and  F in  Lagrange’s  equations 


-'3 


. a.'.- +«?_<■ 

•H\dqJ  dq „ 


the  differential  equations  of  motion  for  each  normal  coordinate  are: 


<?»  + 


n7n-  p . 


4 P 


n-ir-  G 2 „ 

P »• + If  ' > *■  - pT( y" 


(5) 


where  (n  =»  1,  3,  5,  . . .). 

The  solutions  of  these  independent  ordinary  differential  equations,  satisfying  initial 
conditions  in  Eq.  3,  may  be  easily  obtained  by  use  of  Laplace  Transforms  (10).  Then, 
the  displacement  function  for  the  duration  of  the  blast  pulse  is: 

u>(x,f)  = 2 T"  sin  f # sin  £„(f  - r)  dr,  (6) 

n-u.s  <->  21  Jo 
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where  Hk(t)  - -jL/wj*  p(x)  sin  dx. 


. n1  t*  m 

5"  " 8 p ' P 


VpT-5| 

If  the  blast  pulse  is  of  duration  t,  and  w{x,ti),  w(i,tx)  are  the  displacements  and 


2 00 
/ 


«>(x,o  *■  “ 2 


n—  1.3. ft 


velocities,  respectively,  at  t - then  the  solution  for  t > (,  is: 

«-*.)  ^ | (sin  cos  £,^0  + C08^’><l  • \j/\  sin  £„/ 

+ i ^cos  £„*,  - ^ sin  £„/i^  <p  - cos  £,.<J  sin  *—■>  (7) 

where  <t>  - / u’(x,f,)  sin  ~ ,X  dx 

Jo 

^ “So  ' °('X,t S'n  ^ dx' 

Equation  6 has  the  form  of  a transient  forced  vibration  with  damping  while  Eq.  7 
represents  a damped  free  vibration.  Owing  to  the  presence  of  the  elastic  restoring  force 
the  motion  of  the  structure  following  the  pulse  will  consist  of  a damped  oscillation  about 
the  original  configuration  of  equilibrium. 

It  is  interesting  to  examine  two  special  limiting  cases  of  the  preceding  solution. 
The  first  case  is  that  of  a structure  for  which  the  internal  viscous  resistance  is  so  small 
that  it  may  be  neglected;  i.e.,  effectively  a Hooke  solid.  The  solution  then  becomes: 


(r)  sin  pn(t  - t)  dr  for  0 £ l £ 


(8) 


nTX 


, ,,  "V  1 nirx  C u 

u>(x,t)  - Z — sin  --nr  I H, 

u>(x,<)  =*  ^ 2 s'n  P**>  + p cos  P«b  j s'n  P»*  + (^t>  cos  P»b  - sin  pKt\ j cos  p„f  J si 

for  t £ (9) 

Although  this  solution  represents  an  undamped  vibration  of  infinite  duration,  the 
counterpart  of  which  is  not  found  in  physical  structures,  it  may  produce  valuable  informa- 
tion, particularly  in  the  cases  of  brittle  structures  or  structures  with  a definite  yield  point. 
For  such  structures  the  critical  displacement  or  stress  is  known.  The  solution  will  reveal 
whether  the  critical  value  is  reached,  and  if  so,  the  time  at  which  this  occurs,  as  well  as 
fully  describing  the  motion  prior  to  fracture  or  yielding. 

If  the  solution  is  further  particularized  to  the  case  of  the  specific  blast  pulse 


P{x,t)  " PoC~ 


Eq.  8 reduces  to 

w(x,t) 


2 Pod 


n-  1,3,6 


PA  lpn{a‘  + pi) 


- cos  pj  sm  pHt  > sm 


nirx 

21  ' 


Although  theoretically  this  pulse  is  of  infinite  duration,  its  effect  on  the  structure 
soon  becomes  negligible.  For  normal  values  of  a the  exponential  term  rapidly  becomes 
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negligible  as  t increases,  leaving  what  is  sensibly  a free  vibration  of  determinate  amplitude. 

The  other  limiting  case  is  that  of  a structure  of  which  the  elastic  restoring  forces  are 
comparatively  so  weak  that  they  may  be  neglected.  This  amounts  to  setting  the  shear 
modulus  G equal  to  zero,  in  which  event  the  structural  resistance  law  Decomes  that  of  the 
Newtonian  fluid,  Fig.  1(b).  For  this  case  the  general  Eq.  6 and  7 reduce  to: 

= 2 A sin  r*.Mn  " dr  for  0 £ t $ l,  (10) 

— U.s2*-  21  Jo 

w(x,t)  *=  w(x,h)  + y £ 1 1 - sin  for  t £ l,.  (11) 

n-1.3,s""‘  2' 


As  an  illustration  of  the  nature  of  these  solutions,  consider  the  blast  loading  of  a 
structure  having  the  dimensions  75  ft  X 75  ft  X 100  ft  (height)  which  is  located  at  a distance 
of  3600  ft  from  the  ground  zero  of  a “nominal”  atomic  bomb.  Then,  by  approximate 
methods  (3),  the  impulsive  loading  function  may  In?  estimated  to  be: 


PM  =»  P(t) 


0 for  t < 0 

- 324000  M - V50'  for  0 £ t % 0.7 

( 0 for  t > 0.7  sec., 


where  the  origin  for  time  is  taken  as  the  instant  when  the  shock  wave  reaches  the  structure. 
A plot  of  this  loading  function  is  included  in  Fig  4.  The  remaining  necessary  data  were 
assumed  to  be:  .4  = KX)  sq  ft,  weight  = 322  lb  cu  ft,  and  P(t)  = 405,000  lb  sec  sq  ft. 

During  the  non-zero  portion  of  the  pulse  the  displacement  function,  obtained  by 
substituting  these  data  in  Eq.  10,  is 


X | 

tc(z,t)  - 0.14733  ^ sin 

n-l.3,s 


niri  rWnVn1  - 2 )(*-«* 

2 f L 


- n*(13n*  - 24)e-**  + 4(7n« 
(«’  — 2)1 


The  displacements  of  points  originally  on  a vertical  line,  computed  from  this  solu- 
tion, are  shown  graphically  in  Fig.  4 for  successive  instants  of  time  during  the  pulse.  Also 
shown  is  the  velocity-time  variation  of  representative  points  on  the  structure.  It  may  be 
seen  that  the  structure  is  substantially  at  rest  by  the  end  of  the  pulse,  making  unnecessary 
the  determination  of  a separate  solution  for  t > This  feature  is  a result  of  the  particular 
parameters  chosen  for  this  example,  and  generally  would  not  be  the  case.  The  limiting 
displacement,  which  for  all  practical  purposes  is  reached  during  the  pulse,  is: 

w = 1.916  2 ~i  8'n  ~£r 

,-u.8n’  21 


Although  theoretically  the  velocities  do  not  become  identically  zero  in  finite  time 
for  Eqs.  10  and  11,  the  presence  of  internal  friction  in  real  materials  would  suffice  to  bring 
the  structure  to  rest  when  the  velocities  become  sufficiently  small. 

Attention  should  be  called  to  the  basic  difference  between  this  solution  and  that 
for  the  preceding  cases.  In  the  absence  of  elastic  restoring  forces  there  is  no  vibration,  and 
permanent  deformation  of  magnitude  dependent  upon  the  applied  blast  pulse  occurs. 

In  the  preceding  analysis  it  was  assumed  that  the  coefficients  G,  p,  p,  A were  con- 
stants throughout  the  structure.  Solutions  may  be  obtained  for  more  complicated  cases 
where  these  coefficients  vary  with  the  longitudinal  dimension  r by  use  of  numerical  methods. 
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HORIZONTAL  DISPLACEMENT  OF 
VERTICAL  AXIS  OF  BUILDING  (f««») 


When  the  coefficients  are  functions  of  the  displacement  or  its  derivatives  one  encounters  an 
important  class  of  non-linear  problems,  which  are  beyond  the  scope  of  this  study.  How- 
ever, a treatment  will  now  be  given  of  an  important  technological  case  in  which  the  viscosity 
is  distributed  throughout  the  height  of  the  structure  in  a discrete  or  discontinuous  manner. 
In  this  case,  the  structure  is  also  assumed  to  be  rigid  except  at  the  interfaces  of  adjoining 
rigid  blocks. 
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BLAST  LOADING  OF  CEMENTED  RIGID  BLOCKS 


Another  structural  model  of  considerable  interest  is  that  of  a vertical  system  of 
rigid  masses  A/  acted  upon  bv  external  impulsive  forces  P and  coupling  forces  F which  are 
internal  to  the  structure,  as  shown  in  Fig.  5.  In  this  model  it  is  assumed  that  the  coupling 
forces  are  transmitted  by  intervening  material  of  negligible  mass.  This  model  is  especially 
applicable  to  structures  composed  of  massive  blocks  of  masonry  cemented  together;  it 
also  has  application  to  buildings  for  which  the  mass  may  reasonably  lie  assumed  to  be  con- 
centrated at  floor  levels.  The  latter  application  was  suggested  in  Appendix  A of  Ref.  3 
and  was  used  earlier  by  Riot  (11)  for  a building  considered  as  a two-degree-of-freedom  system 
subjected  to  earthquake  motions. 

As  a result  of  the  representation  of  the  structure  as  a system  with  a finite  number  of 
degrees-of-freedom  the  equations  of  motion  are  n ordinary  differential  equations  of  the 
form 

MrW,  - Fr  ,4-  Fr  = Pr(t),  (12) 

where  r = 1,2,3,  . . .,  n 

In  the  notation  employed  (Fig  5),  the  subscript  letters  identify  a particular  mass 
and  the  coupling  material  immediately  below  it  The  displacement  functions  u>i(t), 
u\(l),  . . .,  Un(t)  are  used  to  define  the  motions  of  the  n masses.  The  sign  convention  for 
positive  forces  and  displacements  is  that  indicated  in  Fig.  5.  The  boundary  conditions  for 
this  structure  may  be  expressed  as: 

F0  = 0,  = 0. 


The  coupling  forces  F may  be  of  many  types,  some  of  the  most  important  being 
those  based  on  the  resistance  iaws  depicted  in  Fig.  1.  In  the  case  where  the  coupling 
forces  result  from  the  elastic  shear  deformation  of  the  coupling  material,  Kq.  12  becomes: 


II  *'I|  -\Gr  — I , ( df-lfif-l  3 r (l  r i, 

.17, U', ^ U\-l  + I £ + -jj—  1 U- , UV-fl  = Pr(t), 


(13) 


where  r = 1,2,  . . .,  n 

A = effective  shearing  area  of  designated  coupling  material 
G = effective  shear  modulus  of  designated  coupling  material 
h = height  of  designated  coupling  materia! 

The  solution  of  Kq.  13  may  Ik*  obtained  by  any  of  the  standard  methods,  such  as 
use  of  the  operational  calculus,  at  least  for  a specific  problem. 

This  solution  will  have  the  form  of  an  undamped  vibration  of  each  element  of  the 
system.  The  remarks  made  pertaining  to  the  value  of  the  solution  for  the  undamped 
continuous  structure  apply  equally  well  to  this  solution. 

The  case  of  the  discrete  system  with  coupling  forces  of  the  St.  Venant  body  type 
(Coulomb  friction  or  “perfectly  plastic  material”)  has  been  considered  in  a previous 
paper  (2)  and  also  examined  with  greater  thoroughness  in  an  unpublished  essay  (12). 
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Coupling  forces  of  the  viscous  type  may  also  occur,  as  exemplified  by  the  case  of  the 
structure  composed  of  massive  blocks  cemented  together  by  an  essentially  viscous  material. 
For  such  structures  Eq.  12  becomes: 


1#  Hr—\Ar—\  , / M'-I^  '-I  l _ D f,\  /iA\ 

M,Wr u>,_|  + ( J to, UWl  “ r r(l),  (14) 


where  r - 1,2,3,  . . .,  n. 

Solutions  of  Eq.  14  may  be  easily  obtained  for  small  values  of  n;  using  the  notation 


C, - 


Hr  A, 

K ’ 


where  r - 1,2,3,  . . .,  n,  the  solution  for  the  single-degree-of-freedom  system  is: 
and  for  the  two-degree-of-freedom  system  is: 


— JW*i(A 


(O  - 6)(C,  ± Ct) 
C,C, 


)■ 


,<«+4)<l-r) 


1 / 1 + (o  + 6)(C, 


-±(± 

2 b\Mx 


o * / / 

+ 2^r  )26+  (a  - b)eim**u,~T)  - (a  + 6)e(,-*l<'-'' 1 />»(r)Jdr  + tci(lo) 

"'(%r,)'i,'(W  + £HI'  - K1  - 1)'“**'""“'  - 


05) 


,(•—4)  (I— to)  1 


, u>l(k)  I (,+*)(!_„)  _ (.-*>(,-411 

26  1 ' 


(16) 


u>i(f)  - ~^/'<[j26+  (a  - 6)eu+kl“-')  - (a  + 6)c(-*>(‘-”|/>.(t) 

+ (26  + + a - b)r“+*"-''  - + a + 6^ «<•-*>"-'' j P,(t)J  dr 

+ write)  + {^’  «0,«.)  + £ ii>,(fo)}{  1 - |(l  - ?)«*•*»-» 


where  a 


(3/,  + M\)C\  + MxCr 
2M  |A/| 


(17) 


V^Af,  + A/,)*C*  + 2Ml(Ml  - Mt)CxCr  + 3/?Cj 

2M\M\ 

It  frequently  happens  that  the  external  load  function  has  different  analytical  defini- 
tions for  succeeding  intervals  of  time  (Fig.  6). 

In  such  cases  it  is  necessary  to  obtain  solutions  for  each  intervul,  using  the  final 
conditions  of  displacement  and  velocity  for  one  interval  as  initial  conditions  for  the  next 
interval.  The  solutions  for  one  and  two  degives-of-freedom  an?  presented  in  a form 
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where  initial  conditions  are  specified  at  1 ■»  <o,  an  arbitrary  time,  making  possible  the  direct 
determination  of  the  solution  for  each  interval  with  respect  to  the  same  origin  for  time. 
Obvious  simplifications  can  be  made  when  the  initial  conditions  are  specified  at  1 m 0 or 
when  the  initial  displacement  or  velocity  is  sero.  Equations  15,  16,  and  17  also  give  the 
solution  for  the  period  after  the  impulse  has  ceased,  provided  the  upper  limit  of  integration 
is  taken  to  be  t,,  where  li  is  the  time  of  the  end  of  the  impulse. 

As  an  illustrat^n  of  the  nature  of  the  response  of  a two-degree-of-freedom  system 
with  coupling  forces  of  viscous  type,  Fig.  7,  consider  as  an  example  a system  having  the 
simple  parameters: 

A/,  - 3 C,  = 1 
Mt  = 5 Ci  = 2, 


which  is  subjected  to  the  triangular  blast  pulses  Pi(<),  Pj(0  shown  in  Fig.  8.  These  pulses 
may  be  described  analytically  by  the  functions: 


f 0 for  1 < 0 

P,(0  - j 10  - 0.51  for  0 £ t £ 20 
[ 0 for  1 > 20, 


/MO  = 


0 

24  — 31 
0 


for  1 < 0 
for  0 S 1 £ 8 
for  t > 8. 


Introducing  these  data  into  Eq.  16  and  17,  the  complete  displacement-time  relations 
are  found  to  be: 


tOj(0  “0  1 0 

= 252.56c"  17‘"  - 3.919c-  - 1.1261’  +-  41.511  - 248.65  0 3 1 3 8 

- - 5.29c  + 330.50c-  7‘7*'  - 0.376!7  + 19.001  - 54.39  8 £ 1 3 20 

- - 756.14c"  ,7*“  - 367329c-  7M,‘  + 197  73  1 £ 20 

te,(l>  - 0 15  0 

- 1 19.26c-  ,7*"  +-  4.98c-  7‘7*  - 0.875 17  4-  24.76!  - 124.24  0 £ 1 £ 8 

- - 2.529c-  17»"  - 424.52c- 7‘7«  - 0.125!’  +-  0.75 1 + 2.60  8 3 1 £ 20 

- - 356.75c  - 17‘"  +-  486494c-  7»”'  +-  97.92  < £ 20 


2 

UJ 

U 

< 


Ol 

<✓> 


Fig.  9 


ORO-T-198 


15 


These  displacement-time  functions,  and  the  velocity-time  functions  which  may  be 
easily  obtained  by  differentiation,  are  presented  graphically  in  Fig.  9.  The  maximum  dis- 
placements shown  in  this  figure  are  approached  asymptotically  as  time  increases  indefinitely. 
They  serve  as  an  upper  bound  which  would  be  reached  eventually  if  there  were  no  other 
resisting  forces  present.  As  might  be  expected,  this  solution  has  much  in  common  with 
that  for  the  continuous  structure  having  internal  viscous  resistance  alone,  Fig.  4.  It  is 
also  interesting  to  compare  the  results  of  Fig.  9 with  those  presented  in  “ Blast  Loading  of 
Structures”  (2),  in  which  an  example  with  the  same  parameters  was  solved  for  a system 
having  the  resistance  law  of  the  modified  St.  Venant  body,  Fig.  1(c). 
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DISCUSSION  AND  CONCLUSIONS 


The  type  of  analysis  presented  in  this  paper  is  considered  to  be  of  value  in  predicting 
damage  causer!  by  blast  against  a multi-storied  building  or  similar  structure.  This  point 
of  view  is  advanced  in  spite  of  obvious  defects  that  exist  in  idealizations  of  the  type  in- 
volved. As  a starting  point  some  such  analysis  is  essential.  Statical  analysis  gives  prac- 
tically no  useful  information  about  dynamic  strains  and  deflections.  Also,  any  analysis 
which  assumes  that  the  structure  is  rigid  during  the  first  portion  of  the  time  of  application 
of  the  load  is  misleading.  This  assumption  is  sometimes  tempting  in  order  to  apply  simpli- 
fied theories  of  plasticity.  An  essential  difficulty  in  such  approach  is  the  lack  of  knowledge 
of  the  distribution  of  momenta  caused  by  the  first  phase  of  the  loading. 

The  Biot  method  (1 1)  of  treating  shear-frame  buildings  may  be  utilized  in  the  analysis 
of  similar  structural  units  occurring  in  equipment  and  machines  subjected  to  dynamic  loads. 
However,  in  order  to  make  the  analysis  more  realistic  it  will  be  necessary  to  determine 
suitable  stress-strain  laws  in  the  particular  cases  of  interest. 

Investigations  involving  the  firmo-viscous  hypothesis  emphasize  the  need  for  more 
reliable  information  on  structural  damping  While  many  tests  have  boon  conducted  to 
determine  the  damping  properties  of  materials  the  use  of  damping  laws  in  structural 
analysis  is  still  in  a rather  embryonic  stage.  Considerable  effort  to  remedy  this  situation 
seems  to  be  warranted  at  the  present  time 
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